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We show that, when graphene is subjected to an appropriate one-dimensional external periodic 
potential, additional branches of massless fermions are generated with nearly the same electron- 
hole crossing energy as that at the original Dirac point of graphene. Because of these new zero- 
energy branches, the Landau levels at charge neutral filling becomes 'i{2N + l)-fold degenerate 
(with A'=0, 1,2, tunable by the potential strength and periodicity) with the corresponding Hall 
conductivity a^y showing a step of size 4(2A^ + 1) e^ /h. These theoretical findings are robust against 
variations in the details of the external potential and provide measurable signatures of the unusual 
electronic structure of graphene superlattices. 



The physical properties of graphene [J, 0, 0] are cur- 
rently among the most actively investigated topics in con- 
densed matter physics. Graphene has the unique feature 
that the low-energy charge carriers are well described by 
the two-dimensional (2D) massless Dirac equation, used 
for massless neutrinos, rather than by the Schrodingcr 
equation [2, ji]. Moreover, graphene is considered to be 
a promising candidate for electronics and spintronics ap- 
plications [4|. 

It has been shown that, because of their gapless en- 
ergy spectrum and chiral nature, the charge carriers in 
graphene are not hindered by a slowly varying electro- 
static potential barrier at normal incidence ^] , analogous 
to the Klein tunneling effect predicted in high-energy 
physics. Direct evidences of Klein tunneling through a 
single barrier in graphene [5j] have been observed in re- 
cent experiments [6|, lZ[ . 

Application of multiple barriers or periodic poten- 
tials, either electrostatic [1, |2, [13, [HI, [I3| or mag- 
netic 0, [a, [11,0, to graphene has been shown to 

modulate its electronic structure in unique ways and lead 
to fascinating new phenomena and possible applications. 
Periodic arrays of corrugations [3, [3] have also been 
proposed as graphene superlattices (GSs). 

Experimentally, different classes of GSs have been fab- 
ricated recently. Patterns with periodicity as small as 
5 nm have been imprinted on graphene thr oug h electron- 
beam induced deposition of adsorbates [20|. Epitax- 
ially grown grap hene on the (0001) surface of ruthe- 
niiuii [2l|j23jS^ and that on the (111) surface of 
iridium 1271 . |28| also show superlattice patterns with 
~3 nm lattice period. The amplitude of the periodic po- 
tential applied to graphene in these surface systems has 
been estimated to be in the range of a few tenths of an 
electron volt [1^. Fabrication of periodically patterned 
gate electrodes is another possible way of making GSs 
with periodicity close to or larger than ~20 nm. 

The quantum Hall plateaus in graphene take on the un- 
usual values of 4(? 4- 1/2) jh where I is a non-negative 
integer [i^. The factor 4 comes from the spin and val- 



ley degeneracies. In bilayer graphene, the quantum Hall 
plateaus are at Ale^/h with I a positive integer [30| . 
These unconventional quantum Hall effects have been ex- 
perimentally verified [3, H, [3l[, providing evidences for 
2D massless particles in graphene and massive particles 
in bilayer graphene. 

In this paper, we investigate the LLs and the quan- 
tum Hall effect in GSs formed by the application of a 
one-dimensional (ID) electrostatic periodic potential and 
show that they exhibit additional unusual properties. We 
find that, for a range of potential shapes and param- 
eters, new branches of massless fermions are generated 
with electron-hole crossing energy the same as that at 
the original Dirac point of pristine graphene. These ad- 
ditional massless fermions affect the LLs qualitatively. 
In particular, the LLs with energy corresponding to the 
Fermi energy at charge neutrality (i. e. , zero carrier den- 
sity) become 4(2iV-h l)-fold degenerate {N = 0, 1, 2, ...), 
depending on the strength and the spatial period of the 
potential (pristine graphene corresponds to ~ 0). 
Accordingly, when sweeping the carrier density from 
electron-like to hole-like, the quantum Hall conductivity 
in such a GS is predicted to show an unconventional step 
size of 4(2A^ -t- 1) /h that may be tuned by adjusting 
the external periodic potential. 

In our study, the electronic structure of the GSs is eval- 
uated using the methods developed in Ref. |3; we evaluate 
the bandstructure of the GS numerically by solving the 
2D massless Dirac equation with the external periodic 
potential included using a planewave basis. Similarly, 
to obtain the LLs, the eigenstates of the GSs under an 
external perpendicular magnetic field are expanded with 
planewaves. We work in a Landau gauge with the vec- 
tor potential depending on the position coordinate along 
the direction of the periodicity of the GS, and a zigzag 
form for the vector potential with a very large artificial 
periodicity (large compared to the GS periodicity) is em- 
ployed to mimic a constant magnetic field near the origin 
in position space (sl] . We have checked that the LLs are 
converged in energy to within less than 1 % with respect 
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to the size of the supercell for the vector potential and 
the kinetic energy cutoff for the planewaves. The size of 
the largest supercell and that of the smallest sampling 
distance in real space used are 400 and 0.05 in units of a 
single unit cell, respectively. 

Figure [T{a) shows a GS formed by a Kronig-Penney 
type of electrostatic potential periodic along the x direc- 
tion, with lattice parameter L and barrier width L/2. Re- 
markably, unlike that in graphene [Fig. [Tfb)], the band- 
structure in a GS [Fig. [TJc)] can have, depending on the 
potential barrier height J7o, more than one Dirac point 
with fca; = having exactly the same electron-hole cross- 
ing energy [s^l. As Fig.[ljc) shows, the number of Dirac 
points for this type of GSs increases by two (without con- 
sidering the spin and valley degrees of freedom) whenever 
the potential amplitude exceeds a value of 



= AnNhvo/L 



(1) 



with N a positive integer. The value of the potential 
barrier given in Eq. ([T]) corresponds to special GSs in 
which the group velocity along the ky direction vanishes 
for charge carriers whose wavevector is near the origi- 
nal Dirac co ne [ e. g. , the Dirac cone at the center in 
Fig. [Ijc)] [H,!!^. All the findings in this study apply in 
general to GSs made from a periodic potential which has 
both even and odd symmetries, like a sinusoidal type of 
potential. The results for GSs whose odd or even sym- 
metry is broken are discussed in Ref. |36|. 

Figure [5] shows the evolution of the energy of the elec- 
tronic states with = for a GS depicted in Fig. [T] for 
several different values of Uq. As stated above, the group 
velocity along the ky direction becomes zero near ky = 
when the barrier height is given by Eq. ^ [Figs, ^c) 
and^e)]. When Uq has a value between those specific 
values, the position of the additional new Dirac points 
move away from the ky = point along the ky direc- 
tion with increasing Uq. The complex behavior of the 
zero-energy Dirac cones revealed by our numerical calcu- 
lations cannot be derived using perturbation theory 
because ky is not small compared to the superlattice re- 
ciprocal lattice spacing 2tt/L. Moreover, the pseudospin 
character of these additional massless fermions [e. g. , the 
left and the right Dirac cones (not the center one) in 
Fig. [ijc)] are different from that of the original massless 
Dirac fermions. For example, backscattering amplitude 
due to a slowly varyingpotential within one of the new 
cones does not vanish [36|. 

A natural question arising from this peculiar behav- 
ior in the electronic structure of a GS, which is topo- 
logically different from that of pristine graphene, is how 
the LLs are distributed. Figure [3] shows the calculated 
LLs of the ID Kronig-Penney GSs depicted in Fig. [T] 
for various values of Uq [s^I- When the superlattice po- 
tential modulation is moderate [Fig. [IJb)] , the spacings 
between neighboring LLs become smaller than those in 
pristine graphene [Fig. [31^ a)], owing to a reduction in the 



band velocity. Once Uq becomes larger than Air hvo /L {= 
0.4 eV for L = 20 nm), the zero-energy LLs (correspond- 
ing to zero carrier density) become three-fold degener- 
ate [Fig. md)]. An important point to note is that this 
degeneracy is insensitive to Uq over a range of Uq near 
6tt Uvq I L because the topology of the electron bands does 
not change with this variation 3^, 4^. Moreover, even 
though the massless particles of the different Dirac cones 
may have different band velocities, the degeneracy of the 
zero-energy LLs is not affected. 

The dependence of the Hall conductivity <Jxy on the 
charge carrier density n most directly reflects the degen- 
eracy of the LLs. Figure d] schematically shows that, 
depending on the superlattice potential parameters, Gxy 
of the GSs considered has a 4(2A'' -I- 1) /h step as the 
density is scanned from hole-like to electron-like carri- 
ers. (We have put in the additional factor 4 coming from 
the spin and valley degeneracies in this discussion and in 
Fig-El) Because the degeneracy of the LLs in the ID GSs 
is insensitive to a variation in Ua, this qualitative differ- 
ence in Uxy of the ID GSs from that of pristine graphene 
(Fig. [4]) is expected to be robust, and will provide a mea- 
surable signature of the unique electronic structure of the 
ID GSs. 

In conclusion, we have shown that the electronic struc- 
ture of ID graphene superlattices can have additional 
Dirac cones at the same energy as the original cones at 
the K and K' points of pristine graphene. These new 
massless particles contribute to a 4(2A^ -I- l)-fold degen- 
eracy in the zero-energy Landau levels, whose signature is 
reflected in a 4(2A^-|-1) /h Hall conductivity step where 
TV = 0, 1, 2, ... depending on the superlattice potential 
parameters. This feature of the electronic structure of 
the ID graphene superlattices gives rise to new proper- 
ties for the quantum Hall effect. Equally importantly, 
these new phenomena may provide a direct way to char- 
acterize the peculiar electronic structure of these systems 
experimentally. 
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Note added. — After submission, we became aware of a 
recent theoretical work l f4]| confirming the newly gener- 
ated massless fermions reported in this manuscript, with 
applications to transport properties. 
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FIG. 1: (color online) (a) Schematic diagram of a Kronig-Penney type of potential applied to graphene with strength f/o/2 
inside the gray regions and —Uo/2 outside with lattice period L and barrier width L/2. (b) Electron energy in units of el 
(= hvo/L; for example, if i = 20 nm, el = 33 meV) versus wavevector near the Dirac point in pristine graphene. (c) The same 
quantity as in (b) for a GS with Uo = 6ttel- (d) Number of Dirac points (not including spin and valley degeneracies) in a GS 
versus Uo- 





1 


(a) ' ' ^ 


1 'Uo = 0'- 


-(d) 




= 671 / 




s 




-1 


























1 


(b) \ 


1 ^/o = 27^- 


(e)' \ 




= 8%^ j 








-1 
















/ 


\ ^ 








1 










= 1071 A 








-1 
















(c) \^ 











-3-2-10123-3-2-10123 
k.. (2jt/L) k {2n/L) 



FIG. 2: (color online) Electron energy (in units of el = 
hvo I L) versus ky with fc^, = in GSs shown in Fig. [1] for 
several different values of barrier height Uo (specified in each 
panel in units of el). 



SUPPLEMENTARY INFORMATION 

1. Sinusoidal superlattice 

In the main manuscript, we state that the essential fea- 
tures in the electronic structure of graphene superlattices 
(GSs) revealed by considering the Kronig-Penney type of 
potential remain valid for GSs made with different types 
of periodic potentials. In this section, we support this 
by showing the results for GSs with sinusoidal types of 
external periodic potentials [Fig. [Sl^a)]. 

The function a{x) defined by Eq. (6) in Ref. [l3 for 
a sinusoidal type of external periodic potential V{x) = 
Vb sin(27ra;/L) is a(x) = —VoL/Trhvo-cos{2Trx/L). There- 
fore, as shown in Eqs. (9) and (15) of Ref. [H, the 
group velocity at the original Dirac point perpendic- 
ular to the periodic direction is given by Vy — /qWoj 
where vq is the group velocity in pristine graphene and 
/o — JoiLVo/Trhvo). Here, Jo{x) is the zeroth order 
Bessel function of the first kind. Our calculations show 
that a new pair of massless Dirac points are generated 
whenever /o — 0, i.e., Vb is equal to 

Vq = TTXo, N — , (2) 

where xq.at is the A^-th root of Jo{x) (e.g., xq^i = 2.405, 
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FIG. 3; (color online) Landau level energy Ei (in units of 
eb = hvo/ls with Ib ~ y^hc/eB) versus the Landau level 
index i {i = 0, ±1, ±2, ... ) in GSs formed with a ID Kronig- 
Penney potential for several different values of barrier height 
Uo, with lattice period L = 0.51b- The LLs now have a fi- 
nite width AE (shown not to scale and exaggerated in the 
figure) arising from the ky dependence of the energy of the 
electronic states in a perpendicular magnetic field ^3]. Note 
the 3-fold and the 5-fold degeneracies around -Ei = in (d) 
and (f), respectively. (If the spin and valley degeneracies are 
considered, those become 12-fold and 20-fold, respectively.) 



xo,2 = 5.520, etc.) [Fig.l^b)]. 

Figure [5fc) shows the energy bandstructure of a sinu- 
soidal type of GS with Vb = 4.07r • ^ . Because this value 
of Vq is between Vq and Vq, a pair of new zero-energy 
massless Dirac cones are generated, and they clearly af- 
fect the Landau level degeneracy [Fig. [DJd)] in the same 
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FIG. 4: (color online) Hall conductivity o-^y versus carrier 
density (with an artificial broadening for illustration) for a 
ID Kronig-Penney GS with [/q near 67r?ivo /i (solid red line) 
is compared to that of pristine graphene (dashed blue line). 



way as discussed in the main manuscript for a Kronig- 
Penney type of GS. 



2. Effects of symmetry breaking on the 
newly generated massless fermions 



In this section, we discuss the effect of symmetry break- 
ing of the external periodic potential on the newly gener- 
ated massless fermions. (The case of random perturba- 
tion is discussed in Ref. 42 of the main manuscript. In 
this section, we focus on a periodic perturbing potential 
that breaks the even or odd symmetry.) 

Figure [S{a)-(c) repeats the results shown in the main 
manuscript for a Kronig-Penney type of periodic poten- 
tial having both even and odd symmetries. If the odd 
symmetry is broken by adding an appropriate pertur- 
bation [Fig.[S{d)], new branches of massless fermions are 
still generated [Fig.lSfe)]. The energy at these new mass- 
less Dirac points however is different from that of the 
original Dirac point [Fig. [HJe)]. Even though Ei — Q 
Landau level does not have the degeneracy coming from 
multiple Dirac points, some lower- index Landau levels 
still show this kind of degeneracy [Fig. [6l^f)]. 

If the even symmetry is broken through a perturbing 
potential [Fig. [6l^g)], new Dirac points are not generated 
[Fig. EI^li)]. However, the signature of newly generated 
states may still be probed with photoemission experi- 
ments or transport measurements (4l| . 
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FIG. 5: (color online) (a) Schematic diagram of a sinusoidal type of potential applied to graphene with lattice period L and 
potential amplitude Vb [^(a;) = Vo sin(27ra;/L)]. (b) Number of Dirac points (not including the spin and valley degeneracies) in 
a GS versus Vb in units of £l (= hvo/L; for example, if L = 20 nm, el = 33 meV). (c) Electron energy versus wavevector near 
the original Dirac point {k^ = fcy = 0) for a GS with Vb ~ 4.07re_L. (d) Landau level energy Ei (in units of eb = Tivq/Ib with 
Ib = ^Jfic/eB) versus the Landau level index i (i = 0, ±1, ±2, ... ) in a GS formed with a sinusoidal potential with Vb = 4.07r£L 
and L = 0.5/fl. Note the 3-fold degeneracy (becoming 12- fold degeneracy when the spin and valley degeneracies are considered) 
around Ei = Q. 
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FIG. 6: (color online) (a) Kronnig-Penney type of potential V{x) given by f/o/2 for < a; < L/2 and ~Uo/2 for L/2 < x < L 
with lattice period L. (b)Electron energy (in units of ei, = hvo/L) versus ky with fca; = in a GS formed by the periodic 
potential in (a) with Uo — Gtyel (c) Landau level energy Ei (in units of £b = hvo/ls with Ib = y/hc/eB) versus the Landau 
level index i {i — 0, ±1, ±2, ... ) in a GS depicted in (b), with lattice period L — 0.51b- (d) to (f): Same quantities as in (a) to 
(c) for a periodic potential V{x) with a perturbation that breaks the odd symmetry. The perturbing potential AV{x) within 
one unit cell is given by -|-10 % of the potential amplitude (Uo/2) for L/S < a; < 3L/8 and zero otherwise. Blue dashed circle in 
(f) shows a three-fold degenerate set of Landau levels, (g) to (i): Same quantities as in (a) to (c) for a periodic potential V{x) 
with a perturbation that breaks the even symmetry. The perturbing potential /S.V{x) within one unit cell is given by -j-lO % 
and —10 % of the potential amplitude {Uo/2) for L/A < x < L/2 and for L/2 < x < 3I//4, respectively, and zero otherwise. 
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FIG. 7; (a) and (b): Calculated overlap of two quasiparticle 



states lA^kCr) and liP^.yiv), 



i(k'-k).r 



GS depicted in Fig. [S{a) with Uq — 6neL versus 9]^ and 9^/ 
which are the angles between the kx axis and wavevectors k 
and k', measured from the newly-generated massless Dirac 
point (appearing when Uo = 4neL and moving along the ky 
direction as Uo is increased further), respectively. The overlap 
is shown in a gray scale (0 in black and 1 in white). The 
results show negligible dependence on |k| and |k'| when they 
are smaller than ^ 0.05 x 2tt/L (in the figures, |k| = |k'| — 
0.02 x2ir/L). The two states are in the same band (s' = s) in 
(a) and are in different bands (s' — —s) in (b). (c) and (d), 
and, (e) and (f): Same quantities as in (a) and (b) for GSs 
with Uo = Sttel and Uo = IOttel, respectively. 



3. Pseudospins of new massless fermions 

The pseudospin character of the newly-generated 
massless states are different from that of the origi- 
nal Dirac fermions. In order to illustrate the pseu- 
dospin character of these states, we numerically calculate 



the overlap 



of two quasiparticle 



states V's,k(r) and tps',k'{^) in a GS having wavevectors 
k and k' measured from the newly-generated Dirac point 
(appearing when Uq = 47reL and moving along the ky di- 
rection as Uo is increased further). The behavior shown 
in Fig. [3 which corresponds to the overlap of the pseu- 



dospin part of the wavefunctions, is robust if the magni- 
tudes of k and k' are smaller than ~ 0.05 x 2-k/L. 

As mentioned in the main manuscript, the pseudospin 
character of these additional massless fermions [e. g. , the 
left and the right Dirac cones (not the center one) in 
Fig. [He)] are different from that of the original massless 
Dirac fermions. Backscattering amplitude due to a slowly 
varying potential within one of the new cones does not 
vanish (Fig. [7]) . 
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